Exact results for the Barabasi model of human dynamics 
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Human activity patterns display a bursty dynamics, with interevent times following a heavy tailed 
distribution. This behavior has been recently shown to be rooted in the fact that humans assign their 
active tasks different priorities, a process that can be modeled as a priority queueing system [A.-L. 
Barabasi, Nature 435, 207 (2005)]. In this work we obtain exact results for the Barabasi model with 
two tasks, calculating the priority and waiting time distribution of active tasks. We demonstrate 
that the model has a singular behavior in the extremal dynamics limit, when the highest priority 
task is selected first. We find that independently of the selection protocol, the average waiting time 
is smaller or equal to the number of active tasks, and discuss the asymptotic behavior of the waiting 
time distribution. These results have important implications for understanding complex systems 
with extremal dynamics. 

PACS numbers: 89.75.Da,02.50.-r 



Several problems of practical interest require us to un- 
derstand human activity patterns Typical exam- 
ples are the design of telephone systems or web servers, 
where it is critical to know how many users would use the 
service simultaneously. The traditional approach to char- 
acterize the timing of human activities is based in two 
assumptions: the execution of each task is independent 
from the others and each task is executed at a constant 
rate 0, 0, E3, d ■ A specific task, such as sending emails 
or making phone calls, is then modeled as a Poisson pro- 
cess , characterized by a homogeneous activity pattern. 
More precisely the time interval between two consecutive 
executions of a task follows an exponential distribution. 
An increasing amount of empirical evidence is indicating, 
however, that human activity patterns are rather hetero- 
geneous, with short periods of high activity separated by 

long periods of inactivity [UlilaHEHElEl This 

heterogeneity is characterized by a heavy tail in the dis- 
tribution of the time interval between two consecutive 
executions of the given task [j| 0, 0] . 

In practice the execution of one task is not indepen- 
dent for the others. Humans keep track of a list of active 
tasks from where they decide what to do next, the se- 
lection of one task implying the exclusion of the others. 
This picture lead Barabasi to model the task manage- 
ment by a human as a queueing system, where the hu- 
man plays role of the server Queueing systems 
have already received some attention in the physics liter- 
ature [131 Il4l Il5l Il6| . This interest is motivated by the 
observation of a non-equilibrium phase transition from a 
non-congested phase with a stationary number of active 
tasks to a congested phase where the number of active 
tasks grows in time. In the non-congested phase the mean 
waiting time before the execution of an active task is fi- 
nite. When approaching the phase transition point the 
mean waiting time diverges, while it grows with time in 
the congested phase. 

The Barabasi model belongs, however, to a new class 



of queueing models with a fixed number of active tasks. 
In this case the behavior of interest comes from the task 
selection protocol. In the extremal dynamics limit, when 
the highest priority task is selected first, numerical sim- 
ulations and heuristic arguments show that most of the 
tasks are executed in one step, while the waiting time 
distribution of tasks waiting more than one step exhibits 
a heavy tail [jj. Yet, further research is required to ob- 
tain the scaling behavior in the vicinity of this singular 
point. 

In this work we obtain exact results for the Barabasi 
model, allowing us to prove previous conjectures based 
on heuristic arguments and numerical simulations, and 
creating a solid background for future research. We cal- 
culate the priority and waiting time distribution of those 
tasks remaining in the list for the case of two active tasks. 
We corroborate the observation of a singular behavior in 
the limit when the task with the highest priority is se- 
lected first, and derive the corresponding scaling behav- 
ior. We also obtain an upper bound for the average wait- 
ing time, which is independent of the selection protocol. 
Based on this result we discuss the asymptotic behaviors 
of the waiting time distribution. All the results presented 
here were checked by numerical simulations, providing a 
perfect match with the theoretical curves. 

Barabasi model: The Barabasi model is defined as fol- 
lows. A human keeps track of a list with L active tasks 
that he/she must do. A priority x > is assigned to 
each active task when it is added to the list, with a prob- 
ability density function (pdf) p{x). The list is started at 
t = by adding L new tasks two it. At each discrete time 
step t > the task in the list with the highest priority 
is selected with probability p, and with probability 1 — p 
a task is selected at random. The selected task is exe- 
cuted, removed from the list, and a new task is added. 
The control parameter p interpolates between the ran- 
dom selection protocol at p — and the highest priority 
first selection protocol at p = 1 . 
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The numerical simulations indicate that the case L = 2 
already exhibits the relevant features of the model 
Furthermore, if we focus on a single task, such as send- 
ing emails, we can model the active tasks list as a list 
with two tasks, one corresponding to sending emails and 
the other to doing something else. Within this scenario 
the waiting time coincides with the time between two 
consecutive executions of the corresponding task. Thus, 
the L — 2 case provides us with a minimal model to 
study the statistical properties of the time between the 
consecutive execution of specific tasks. 

Consider the Barabasi model with L = 2. The task 
that has been just selected and its priority has been 
reassigned will be called the new task, while the other 
task will be called the old task. Let p(x) and R(x) = 
Jq dxp(x) be the priority pdf and distribution function 
of the new task, which are given. In turn, let pi(x, t) and 
Ri(x,t) = J* dxpi(x, t) be the priority pdf and distri- 
bution function of the old task at the i-th step. At the 
(t + l)-th step, there are two tasks on the list, their pri- 
orities being distributed according to R(x) and Ri(x,t), 
respectively. After selecting one task the old task will 
have the distribution function 




FIG. 1: Old task priority distribution for the case of a uniform 
new task priority distribution function, p(x) = 1 and R(x) = 
x in < x < 1, as obtained from (lines) and numerical 
simulations (points). The case p = corresponds with the 
random selection protocol with Ri(x) — R(x) — x. 
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Rl(x,t+1)= / dx'p 1 (x',t)q(x')+ dx'p(x)q 1 (x',t) , 
Jo Jo 

where 



q(x) = p[l-R(x)] + (l-p)- 



(2) 



is the probability that the new task is selected given the 
old task has priority x, and 



q 1 (x,t)=p[l-R 1 (x,t)] + {l-p)^ (3) 

is the probability that the old task is selected given the 
new task has priority x. In the stationary state, Ri (x, t+ 
1) = Ri(x,t), from we obtain 



i+p 

2p 



(4) 



To analyze this result let us consider its limiting cases. 
When p — > J3J results in 



i.e. pi(x) is concentrated around x — 0. This result 
implies that in the limit p — > 1 the new task will always be 
selected for execution, resulting in a waiting time r = 1 . 
We are going to return to this result after computing the 
waiting time distribution. The progression between these 
two limiting cases is illustrated in Fig. where we plot 
i?i (x) Q as a function of x for a uniform distribution in 
< x < 1. As p increases from zero Ri(x) deviates more 
from R(x), resulting in a higher probability to obtain 
smaller values of x. When p — 0.999, R\{x) grows to a 
value of almost one in a very short x-range close to x — 0, 
approaching the limit distribution © • 

Next we turn our attention to the waiting time dis- 
tribution. Consider a task with priority x that has just 
been added to the queue. The selection of this task is 
independent from one step two the other. Therefore, the 
probability that it waits t steps is given by the product 
of the probability that it is not selected in the first r — 1 
steps and that it is selected in the r-th step. The prob- 
ability that it is not selected in the first step is q\{x), 
while the probability that it is not selected in the subse- 
quent steps is q(x). The integration over the new task's 
priorities results in 



lim Ri(x) 

p^O 



R(x) 



(5) 



Indeed, this limit corresponds to the random selection 
protocol and, therefore, the priority distribution of old 
tasks is equal to that of new tasks. On the other hand, 
when p — > 1 from W e obtain 



Pr 



J™dR(x)[l- qi (x)} , T = \ 

/ °° dR(x) qi (x) [1 - q{x)\ q( X y- 2 , r > 1 



Using ancl integrating Q we finally obtain 



(7) 
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FIG. 2: Waiting time probability distribution for the case 
of a uniform new task priority distribution, p(x) = 1 and 
R(x) = x in < x < 1, as obtained from iJHJ (pluses) and 
numerical simulations (open symbols). The inset shows the 
fraction of tasks with waiting time r = 1 as a function of p, as 
obtained from JHJ (line) and numerical simulations (points). 



FIG. 3: Average waiting time of executed tasks vs the list 
size for the case of a uniform new task priority distribution 
function, p(x) = 1 and R(x) = xinO<x<l, as obtained 
from 1161 (lines) and numerical simulations (points). 



where 



Pr = 



l _ i_£l ln l±£ 

4p l—p 



4p 



m -( 



T = 1 



r > 1 



(8) 

Note that P T is independent of the p(x). Indeed, what 
matters for a task selection is its relative priority with 
respect to other tasks, resulting that all dependences with 
x in ©-(0J and J7J) appears via R(x). 

As before, let us consider the limiting cases. In the 
limit p — > from JBJ it follows that 



lim P T = 



(9) 



for r > 1. This limit corresponds with the random selec- 
tion protocol, where a task is selected with probability 
1/2 on each step. In the other limit, p — > 1, we obtain 



l + 0(i^ln(l-p)) , r = l 



lim P T = 



(10) 



r > 1 



In this case almost all tasks have a waiting time r = 1, 
while the waiting time of tasks that are not selected in 
the first step follows a power law probability distribution. 
This picture is illustrated by a direct plot of P T in JHJl. 
In Fig. we plot P T vs r for a uniform distribution in 
< x < 1. For p — 0.9 the probability distribution Pi- 
is dominated by an exponential cutoff. This exponential 
cutoff can be derived from (|SJ) by taking the limit t — > oo 
with p fixed, resulting in 



Pr ~ ; cxp 

At V t 



(11) 



r = In 



1+p 



(12) 



When p — > 1 we obtain that To — > oo and, therefore, 
the exponential cutoff is shifted to higher r values, while 
the power law behavior P T ~ 1/r becomes more evident. 
The P T vs r curve systematically shifts, however, to lower 
values for r > 1, indicating that this power law applies to 
a vanishing task fraction (see Fig. [2] and (|lll) ). In turn, 
Pi — > 1 when p — > 1, as it is corroborated by the direct 
plot of Pi as a function of p (see inset of Fig. |2J). 

Another characteristic magnitude of a queueing sys- 
tem is the average waiting time of an active task before 
its execution. For L = 2 we can calculate the average 
waiting time directly from ©, obtaining 



2 , 0<p< 1 
1 , P= 1 • 



(13) 



This average is restricted to those tasks that are executed 
and, therefore, for p — 1 we are excluding the task that 
remains indefinitely in the queue. As we show next, we 
can extend this result for L > 2 using simple arguments. 

On each step there are L task in the queue and one of 
them is executed. Therefore 



L-l 



= Lt 



(14) 



where Tj is the waiting time of the task executed at the 



i-th step and tL 



1, . . . , L — 1, are the resident times 



of the L—l tasks that are still active at the t step. From 
((Tljl it follows that 
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1 ' 1 L_1 
(t) = lim — > Tj = L — lim — > t' 



(15) 



For < p < 1 the numerical simulations indicate that all 
active tasks are expected to be executed Therefore 
(t 1 ) < (t) and the last term in (|15fl vanishes when t — > 
oo. In contrast, for p = 1 the numerical simulations 
indicate that after some transient time the most recently 
added task is always executed, while L — 1 tasks remain 
indefinitely in the queue. In this case r[ ~ t when t — > oo 
and the last term in (|15J) is of the order of L — 1 when 
t — ► oo. Based on these arguments we conjecture that 
the average waiting time of executed tasks is given by 



(r) = 



L , <p< 1 
1 , p= 1 . 



(16) 



This results was proved for L = 2 and for L > 

2 it is corroborated by the numerical simulations (see 
Fig. |2| • It is important to note that the equality in 
(|15fl is independent of the selection protocol, allowing us 
to reach more general conclusions beyond the Barabasi 
model. Using (|15|) we obtain 



(r)<L 



(17) 



From this constraint it follows that P T must decay faster 



than 



when t — ► oo. Thus, when r — > oo either 



P T ~ CLT- a , 

where a is a constant and a > 2, or 



Pr 



a f(- 



\ T 



(18) 



(19) 



where tq > and f(x) = 0{bx a ~ 2 ) when x — > oo, where 
b is a constant. For instance, for the Barabasi model 
with L = 2 and < p < 1, a = 1 and /(x) decays 
exponentially (|llfl . in agreement with l|19|l . 

The empirical evidence 0, Hoi Hl| is in favor of the 
second scenario with a = 1. This observation is in agree- 
ment with our expectation of a natural cutoff. For in- 
stance, we might go on a trip and not check emails for 
several days, but sooner or later we are going to do it 
because we receive and transmit important information 
using this communication media. This cutoff time is ex- 
pected to be more related to the necessity of performing 
a given task rather than to the interaction with other 
tasks. In this sense, the random selection of a task in the 
Barabasi model could be interpreted as a task whose pri- 
ority suddenly increases to the maximum value because 



the need to execute it after being on the queue for a long 
time. This indicates future directions of research, con- 
sidering the case when the priority of old tasks may also 
change with time • 

The singular behavior of the Barabasi model is a con- 
sequence of the extremal dynamics rule: the task with 
the highest priority is selected first. Therefore, the con- 
clusions obtained here are also relevant to other complex 
system evolving with extremal dynamics 
this more general context the waiting time represents the 
life time before selection, an important quantity in evo- 
lution models |l8l llif and optimization methods based 
in extremal dynamics |2dl| . Further research is required, 
however, to determine the influence of other factors such 
as correlations among neighbors, which are absent in the 
Barabasi queueing model. 

This work was supported by NSF ITR 0426737, NSF 
ACT/SGER 0441089 awards and the James S. McDon- 
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